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sections. The NLO corrections are found to be non-trivial, resulting in a reduction of the asymmetry. 
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I. INTRODUCTION 



In the last years measurements of the spin asymme- 
tries (N = p, n, d) in longitudinally polarized deep- 
inelastic scattering (DIS) have provided much new infor- 
mation on the spin structure of the nucleoli . Recent 
theoretical leading order (LO) and next-to-leading order 
(NLO) §, g |, |, g g |, | analyses of the data sets 
demonstrate, however, that these are not sufficient to ac- 
curately extract the spin-dependent quark (Aq = —q^) 
and gluon (Ag = g~ — g^) densities of the nucleon. This 
is true in particular for Ag(x,Q 2 ) since it contributes 
to DIS in LO only via the Q 2 -dependence of g\ (or A\) 
which could not yet be accurately studied experimentally. 
As a result of this, it turns out that the a;-shape of Ag 
seems to be hardly constrained at all by the DIS data, 
even though a tendency towards a fairly large positive 
total gluon polarization, f} Ag(x, Q 2 = 10 GeV 2 )<ia: ~ 1, 
was found [g, g g g [£], g. The measurement of Ag 
thus remains one of the most interesting challenges for 
future spin physics experiments. In selecting suitable 
processes for a determination of Ag, it is crucial to pick 
those that, unlike DIS, have a direct gluonic contribution 
already at the lowest order. Here, one thinks in the first 
place of high-pT reactions in nucleon-nucleon collisions, 
which have been tremendously important in the unpolar- 
ized case to constrain the unpolarized gluon density. 

At the moment, the most eagerly awaited experimen- 
tal tool for the 'spin physics' community is the RHIC|[lC|l 
collider at BNL, at which first runs in a proton-proton 
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mode with longitudinally polarized beams are expected 
to be accomplished this year. The center-of-mass energy 
for these pp collisions will be ranging between 200 and 500 
GeV, with luminosities (rising with energy) between 320 
and 800 pb _1 , respectively. One expects about 70% po- 
larization for each beam. Such conditions look extremely 
favorable for studying the spin asymmetries for all kinds 
of high-pT PP processes that are sensitive to the gluon 
density, such as hadron, jet, prompt-photon, or heavy- 
flavor production. Hadrons and jets could be the key 
to Ag: at VS = 200(500) GeV, hadrons(and jets) will 
be extremely copiously produced, and the correspond- 
ing cross-sections will show a strong sensitivity to Ag 
thanks to the dominance of the gg and qg initiated sub- 
processes in some kinematical ranges. While the STAR 
detector permits a clear determination of jets, in the case 
of PHENIX, where the limited coverage in rapidity does 
not allow to observe jets, hadrons (and especially pions) 
can be used as jet surrogates O, . 

In order to make reliable quantitative predictions for a 
high-energy process, it is crucial to determine the NLO 
QCD corrections to the Born approximation. Quite in 
general, the key issue here is to check the perturbative 
stability of the process considered, i.e. to examine the 
extent to which the NLO corrections affect the cross sec- 
tions and (in spin physics) the spin asymmetries relevant 
for experimental measurements. Only if the corrections 
are under control can a process that shows good sensi- 
tivity to, say, Ag at the lowest order be regarded as a 
genuine probe of the polarized gluon distribution and be 
reliably used to extract it from future data. Furthermore, 
the inclusion of extra partons in the NLO perturbative 
calculation also allows to improve the matching between 
the theoretical calculation and the realistic experimental 
conditions. 
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The calculation of the NLO QCD corrections to hadron 
production by polarized hadrons is the purpose of this 
paper. Such a calculation needs the one-loop 2^2 
and tree-level 2 — * 3 polarized amplitudes as input. For- 
tunately, these amplitudes are already known |l3|, pi] . 
Furthermore, several independent methods to calculate 
any infrared-safe quantity in any kind of hard unpo- 
larized collision are at present available in the litera- 
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ture H |T|, [L7 
used in Ref. [IS 



, The formalism of Ref. ]16[ has been 
to construct a Monte Carlo code that 
can calculate any jet infrared-safe observable in hadron 
hadron unpolarized collisions and extended, in Ref. [ jl9| , 
to the polarized case (photo-production of jets has been 
studied in ^(J). In the present paper, we extend the 
method of Refs. Jl6|, [l8| to the case of one single hadron 
inclusive observables. As a result, we will present a cus- 
tomized code, with which it will be possible to calculate 
any infrared-safe quantity corresponding to one-hadron 
production to NLO accuracy, for both polarized and un- 
polarized collisions. It is worth noticing that the code 
presented here can be used to compute the 'resolved' 
contribution to hadron photo-production as well. This 
complements the calculation of the 'direct' component 
presented in Ref. pl| . 

This paper is organized as follows: in section |l| we 
describe the formalis m an d the main ingredients of the 
calculation. In section III we study the perturbative sta- 
bility of the NLO results by looking at the scale depen- 
dence and if-factors. In section [V we present the results 



for the asymmetries at RHIC and study the possibility 
of extracting Ag from pp — > tt data. Section [v| contains 
the conclusions. 



II. FORMALISM AND MAIN INGREDIENTS 



The factorization theorem 1 22 allows to write the cross 
section for the one hadron production in hadronic colli- 
sions as 

da pp^hx jdx x dx 2 dz /f 1 (x u jf 2 (z 2 ,mIi) 

h,h,f 

xda flh ~* fx (xi pi,x 2 p 2 ,Ph/z, fiFi, Vff, Hr) 
xD){z,^ FF ), (1) 

were H\ and H 2 are the colliding particles with momen- 
tum pi and p 2 , respectively, h is the outgoing hadron 
with momentum ph and the sum in Eq. (Q) runs over all 
possible initial and final partonic states. The parton dis- 
tributions f i < are evaluated at the factorization scale 
fiFi, the fragmentation functions at the scale /ijrp and 
the coupling constant, appearing explicitly in the per- 
turbative expansion of the partonic cross section, at the 
renormalization scale /i^. The analogous of Eq.(|l|) for 
the polarized cross section is obtained by replacing the 
parton distributions and the partonic cross section by its 
polarized expressions, Aff 1 ' and dAo"^ 2 ^ x , 
tively. 
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FIG. 1: The polarized gluon (left) and quark densities (right), 
as given by the NLO parametrizations that will be used in this 
paper, at the scale Q 2 — 25 GeV 2 . The patterns for the quark 
densities (= Aq + Aq) are the same as those used for the 
gluon. In the case of DS we show only the results for the + 
scenario, the ones for the — scenario are identical. 



As usual, the (longitudinally polarized) asymmetry is 
defined by the ratio between the polarized and unpolar- 
ized cross sections 



dAa 
da 



(2) 



The method we have used for computing the NLO cor- 
rections and implementing the results in a Monte-Carlo 
like code has been introduced and extensively discussed 
in Refs.Jl^, 18 1, and in Ref. |H. We refer the reader to 



those references for details. It is worth noticing that with 
one single code we can compute unpolarized and polar- 
ized cross sections, since the structure of the corrections 
is the same in both cases. 

As stated in the introduction, there is hardly any ex- 
perimental information on the spin-dependent gluon den- 
sity Ag at present. In contrast, the quark densities are 
far better constrained by the existing data from inclusive 
polarized DIS. This is particularly true for the polarized 
u and d quark densities, which come out rather similar 
in all theoretical analyses performed so far. The spin- 
dependent sea-quark distributions seem less well con- 
strained; however, they are of minor importance for our 
studies. 

In our phenomenological analysis we will try to cover 
as much as possible of the wide range of polarized parton 
densities allowed by the present DIS data: this is espe- 
cially relevant for the gluon density, as we will discuss 
below. For the polarized parton densities, we will use 
the six sets of Ref. f| (DS1±, DS2±, DS3±), obtained 
by constraining the first moment of the polarized gluon 
densities in three different ways. The first moment of the 
gluon of DS3 is the largest and about three times larger 
than for DS1, which has the smallest gluon. The index ± 
refers to the different extremes assumptions done on the 
sea quark densities Au = ±Ad. Furthermore, we will also 
use the 'standard' set of Ref. (from now on, referred 
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to as GRSV), for comparison. All these distributions are 
available at both LO and NLO, the latter corresponding 
to the MS scheme used also in our calculation of the NLO 
partonic cross section. 

In this work, we will mostly concentrate on the phe- 
nomenology of pion production at RHIC with a center-of- 
mass energy of = 200 GeV and transverse momenta 
in the region of 4 GeV < pr < 20 GeV. This implies that 
the polarized parton distributions will mainly be probed 
in the x-range 0.05 < x < 0.3 and at typical scales Q 2 
of the order of 25 GeV 2 . Figure |l| shows the NLO po- 
larized quark and gluon densities of the different sets we 
are going to use. It becomes clear that there is indeed 
a wide range of possible gluon distributions compatible 
with present polarized DIS data. As mentioned earlier, in 
the quark sector, most of the distributions are very sim- 
ilar. The spin-dependent sea-quark densities (not shown 
in Fig. [j]) differ more strongly among the various sets, 
but have a small impact on the cross sections in this kine- 
matical region. In conclusion, since the variations in the 
quark sector are much smaller than the ones for gluons, 
we can expect that any differences between predictions 
for the polarized cross sections (or asymmetries) that are 
found when using different polarized parton density sets, 
are to be attributed to the sensitivity of the observable 
to Ag. 

The size of radiative QCD corrections to a given un- 
polarized hadronic process is often displayed in terms of 
a 'if-factor' which represents the ratio of the NLO over 
LO results. In the calculation of the numerator of K 
one obviously has to use NLO-evolved parton densities. 
As far as the denominator is concerned, a natural def- 
inition requires the use of LO-evolved parton densities. 
In the polarized case, a problem arises for such a defi- 
nition. This is particularly true when the cross section 
is dominated by the gluon distribution: Since the data 
hardly constrain the gluon density, very different results 
for Ag can emerge if the fit is performed at LO or at NLO. 
Therefore, the 'if -factor' for one hadron production, de- 
fined using LO parton densities in the denominator could 
be artificially large or small merely from the fact that the 
gluon is at present so ill-constrained, and not because of 
a 'genuine' effect from higher order partonic corrections. 

In order to avoid this problem, we will follow the def- 
inition introduced to perform the analysis of the results 
for jet production in polarized collisions at NLO |Ts|| , i.e., 
by defining the 'if -factor' as the ratio between the NLO 
and the 'Born' cross section, where the latest corresponds 
to the use of NLO-evolved parton densities (and 2-loop 
expression for a s ) when evaluating the lowest-order par- 
tonic cross sections in the denominator. Nevertheless, it 
is important to remember that the 'if-factor' is not a 
physical quantity and just provides a number to 'quan- 
tify' the effect of the higher order corrections. 

In the unpolarized case, we will use the GRV98 parton 
distributions p3fl . Differences of the order of 10% in the 
cross section are observed when the more recent CTEQ 
p4j or MRST|25| distributions are used. Nevertheless, 



since both DS and GRSV sets of polarized distributions 
used GRV98 as the reference set in the unpolarized case, 
either to impose positivity constraints or to define the 
measured asymmetries, we will restrict the analysis to 
the GRV98 set. 

The other input needed in Eq.(|l|) corresponds to the 
fragmentation functions (FF). Several sets of FF have 
been released in the last years [^6[ |2?j, [29|, as the 
result of QCD analyses of all e + e~ — > h available data. 
As it is well known, e + e _ data can provide information 
the combination D q + Dq in the quark sector and only 
a weak constrain on the gluon fragmentation function, 
which enters at NLO or through the Q 2 evolution. A 
comparison of two different analysis can be found in |27j . 
In particular, large differences are observed between the 
BKK(28) and Kretzer[27| distributions for the gluon in 
the whole range of z and in the quark case at large z (see 
Fig. 7 of (H for details). 

The lack of separation between 'favored' and 'unfa- 
vored' quark fragmentation functions, and the large un- 
certainties on the distributions result on the reduction of 
the predictive power of pQCD calculations in hadronic 
collisions. Besides the large differences than can arise 
in the computed cross section when different sets of FF 
are used, it is not possible to compute observables sep- 
arately for, say, 7r + and 7r~: only cross sections for the 
sum 7r + + 7r~ can be computed with the FF obtained 
from fits to e + e~ data. 

An approach to solve the last problem is attempted in 
Kretzer's FF: a (I — z) suppression for the 'unfavored' 
fragmentation functions, with respect to the 'favored' 
ones, is assumed at the low initial scale //q = 0.4 GeV 2 . 
Within this assumption, a full separation between quark 
flavors is provided, showing a good agreement with the 
available semi-inclusive DIS data. Furthermore, the fit 
by Kretzer relies on the same values of Aqcd used by 
the GRV98, DS and GRSV sets, allowing to perform a 
fully consistent analysis in the case of hadronic collisions. 
Therefore, we will present results for hadron production 
at RHIC usingthe set of fragmentation functions pre- 
sented in Ref . p7[ . 

Finally, it should be noticed that an analytical NLO 
calculation exists for the case of unpolarized hadronic 
collisions J|(| since long time. We have compared nu- 
merically to the code presented in |3(| finding excellent 
agreement. This provides a very important test on our 
calculation and, furthermore, on the structure of the code 
that we will also use to compute the polarized cross sec- 
tion. 



III. PERTURB ATIVE STABILITY 

A reliable error estimate on our NLO results requires 
some knowledge on the size of the uncalculated higher- 
order terms. The best we can do, before higher-order 
terms are computed, is to study the dependence of the 
full NLO results on the renormalization and factorization 
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FIG. 2: Scale dependence of NLO and Born unpolarized cross 
sections. The scales are varied around the default value pt as 
indicated in the text. The LO results have been rescaled by a 
factor 0.1. 



FIG. 3: Scale dependence of NLO and Born polarized cross 
sections. The scales are varied around the default value pt as 
indicated in the text. The LO results have been rescaled by a 
factor 0.1. 



scales. 

Although physical observables are obviously indepen- 
dent of (i, theoretical predictions do have such a depen- 
dence, arising from the truncation of the perturbative 
expansion at a fixed order in the coupling constant a s . 
A large dependence on the scales, therefore, implies a 
large theoretical uncertainty. In order to show how the 
scale dependence is substantially reduced once the next- 
to-leading order corrections are included we will compare 
to the Born result (with a very similar scale dependence 
to the LO one). For some observables it might occur 
that there is a partial cancellation in the scale depen- 
dence when all scales are varied simultaneously, either at 
NLO or at the Born level, and a larger dependence is ob- 
served when the scales are varied independently. There- 
fore, we will look at the scale dependence of the cross 
section when the scales are varied either simultaneously 
or independently, with the others fixed at a default value 
in the last case. 

We begin by presenting the results for the production 
of 7r + + 7r~ in unpolarized collisions at RHIC with y/s = 
200 GeV. The kinematics limits of the PHENIX detector 
< 0.35) have been imposed, but very similar results, 
concerning the perturbative stability, are obtained for the 
case of STAR (with \rj\ < 1). 

In Figure || we show both NLO and Born unpolarized 
cross sections computed at different scales. We take as 
the default scales = hfi = Hff = A* = Pt and com- 
pare to the cases when the renormalization and factor- 
izations scales are varied by a factor of 2 (up and down 
when the scales are varied independently, and up when 
they are varied simultaneously) with respect to the de- 
fault one. As can be observed there is a considerable 



reduction in the scale dependence, showing an improve- 
ment in the perturbative stability for all the scales rele- 
vant to this process. 

Nevertheless, even though of the considerable improve- 
ment, it is worth noticing that the scale dependence is 
still rather large at NLO, of the order of 30% or more, 
and certainly larger than for jet production ]l9| . Actu- 
ally, the scale dependence can be much larger for other 
choices of parton distributions and fragmentation func- 
tions. This issue will be studied in more detail in a future 
publication fl34| . 

We move now to the polarized cross section. To avoid 
the proliferation of curves, we concentrate here on the 
case of the DS2+ set of polarized pdfs. Similar results 
are obtained with other distributions. 

As can be observed in Figure ||, the reduction is the 
scale dependence is even more noticeable than in the 
unpolarized case. A similar situation has been already 
observed for other observables in hadron collisions, like 
jets ML prompt-photon j3l[] and heavy quark produc- 
tion |32|], and should be attributed to the less singular 
behavior of both polarized parton densities and partonic 
cross sections. 

We finish by presenting the X-factors for both polar- 
ized and unpolarized collisions in Fig. ||. As discussed in 
the previous section, if-factors are defined by the ratio of 
the NLO and Born cross sections, with all scales fixed to 
fi = pt- As can be observed, the corrections are larger for 
the unpolarized cross section than for the polarized, re- 
sulting in (about) a 30% reduction of the asymmetry, for 
both PHENIX and STAR. The situation is very similar 
to what has been observed in the case of jet production 



in polarized hadronic collisions |19 
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FIG. 5: Ratios of the unpolarized cross section for different 
combinations of initial (left) and final (right) partonic states. 



FIG. 4: Unpolarized and polarized K— factors for the produc- 
tion of 7T+ + tt" at RHIC with y/s~ = 200 GeV. 



IV. ASYMMETRIES AT RHIC 

In this section we study the phenomenology of pion 
production at RHIC, focusing on the possibility to ex- 
tract the polarized gluon density from such a measure- 
ment. 

It is actually possible to anticipate whether the pro- 
cess is expected to be sensitive to the polarized gluon 
distribution by looking at the relevance of the different 
partonic initial states in the unpolarized case, where the 
parton densities are reasonable well known. 

This is shown on the left-hand side of Figure |[ again 
for 7T+ +TT- production at PHENIX with y/s = 200 GeV. 
All scales have been fixed to fi — pt and R corresponds 
to the ratio between the contribution of the cross section 
by a given partonic initial state and the full one, both 
computed to NLO accuracy. As can be observed, the 
cross section is dominated by the qg initial state at low 
Pt , with an increasing contribution of the qq channel for 
larger transverse momentum. Considering the relevance 
of the gluon distribution in the unpolarized cross section, 
we can expect a large sensitivity on Ag in the polarized 
case. The right-hand side of Figure || corresponds to a 
similar ratio, but in this case, to compare the relevance of 
the contribution of the final state (fragmenting) partons. 
Except at very low px, the quark fragmentation com- 
pletely dominates the cross section. It is worth noticing 
that the situation can change considerably when other 
sets of fragmentation functions are used. 

Before moving to the results for the asymmetries it is 
interesting to find out the range in x and z covered by 
RHIC measurement. This can be studied by looking at 
the dependence of the unpolarized cross section on those 
variables, as shown in Figure [| The results correspond 
to the NLO integrated cross section for pt larger than 
4 GeV and 10 GeV. As can be observed, at yfs = 200 
GeV, the cross section is mostly sensitive to the parton 



distributions in the range 0.05 < x < 0.3 and to the frag- 
mentation functions at large 0.2 < z < 0.8. It should be 
noticed that, unfortunately, this is the region where the 
uncertainties on both fragmentation and (unpolarized) 
parton distributions are considerably large. 

We thus turn to the study of the dependence of the 
asymmetries on the polarized parton densities, as a way 
to determine the sensitivity on the polarized gluon dis- 
tribution. 

We begin with the asymmetry for ir + + ir~ production 
at PHENIX with y/s = 200 GeV. Since the coverage in 
rapidity is rather small, we just concentrate on the px 
distribution. The NLO results, obtained using the sets 
of distributions enumerated in Section 2, are shown in 
Fig.0. 

To estimate the minimum value of the asymmetry ob- 
servable at RHIC, we use the well-known formula 

{ Ah Lh) mm = p2 ^T^' ^ 

where C is the integrated luminosity, P is the polariza- 
tion of the beam, and the factor e < 1 accounts for ex- 
perimental efficiencies; a is the unpolarized cross section 
integrated over a small range in transverse momentum 
(pt bin). The quantity defined in Eq. (||) is plotted as 
'error-bars' in Fig. fr| for e = 1, P — 0.7, a pT-bm size of 
2 GeV and two scenarios for the integrated luminosities 
C = 7 pb _1 and 100 pb _1 . From the figure we see that 
the asymmetry defined in Eq. (|J) is measurable if the 
polarized parton densities are as described by all of the 
sets considered here. Furthermore, it is clear that even 
with the low integrated luminosity expected to be col- 
lected during the first run of RHIC with longitudinally 
polarized beams (~ 7 pb _1 ), it will be possibly to rea- 
sonably estimate the size of Ag from the first px bins. 
Notice that the pattern of the asymmetries computed 
with different polarized distributions follows closely the 
one showed in Fig.jjJ indicating an almost linear depen- 
dence on the gluon density and confirming the relevance 
of the qg partonic initial state also for the polarized cross 
section. 
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FIG. 6: Unpolarized cross section distributions in x = y/xiX2 
and z. 
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FIG. 7: NLO asymmetries for n + +ir~ production at PHENIX 
with y/s = 200 GeV. The 'error bars' estimate the mini- 
mum observable asymmetry assuming integrated luminosities 
of C =7 pb- 1 and 100 pb' 1 . 



Smaller asymmetries are obtained for the y/s = 500 
GeV scenario. Nevertheless, the reduction in the asym- 
metry is compensated by a large increase in the cross 
section, allowing the measurement of the asymmetry up 
to larger values of px- 

Using charged hadrons as jets surrogates allows to in- 
vestigate whether searching for particular hadron in the 
final state provides a larger asymmetry. In this paper 
we concentrate in the cases of 7r + and 7r~, within the 
approximation for the fragmentation functions discussed 
before. 

In Fig.|| we show the NLO asymmetries computed with 
the DS sets for PHENIX with = 200 GeV. The es- 
timate for the minimum observable asymmetry has been 
obtained assuming an integrated luminosity of £ = 100 
pb -1 . As can be observed, the asymmetries are con- 
siderably increased (about a factor of 2) with respect 
to the case of ir + + 7r~, while the minimum observable 
asymmetry only increases by, roughly, a y/2 factor. The 
measurement of ir + production could provide an even 
more powerful tool to extract the polarized gluon dis- 
tributions, after the uncertainties on flavor separation of 
the fragmentation functions are reduced. 

Finally, we move to the case of 7r~ production pre- 
sented in Fig.^. The asymmetries turn out to be much 
smaller than in the other cases, reaching a zero around 
Pt = 10 GeV and becoming negative for larger values of 
the transverse momentum of the pions. Furthermore, the 
sensitivity of the asymmetry on Ag is also considerably 
reduced, becoming very difficult to extract it from such 
a measurement. 

It is possible to understand the differences between the 
asymmetries for n + and tt~ production on simple phys- 



ical basis. The main contribution to the polarized cross 
section, at least regarding the terms depending on the 
gluon density, appears with the following combination of 
parton distributions and fragmentation functions 



Ag(AuDZ + AdD* d ). (4) 



In the kinematical range relevant for RHIC measure- 
ments, Au is found to be positive, while Ad nega- 
tive, with |A<i/Aii| < 1. Furthermore, for tt + the 
'favored' and 'unfavored' fragmentation functions obey 
Du > -^d' wnue the opposite occurs for n~ . In the 
case of the ir + +tt~ production, both fragmentation func- 
tions are the same. Therefore, for tt + the combination 
(Au + AdD^) reaches its maximum, thus increasing 
the dependence of the asymmetry on the polarized gluon 
density, while for 7r~ can be zero, becoming the asymme- 
try dominated by other initial states, like qq. 

Finally, it should be remarked that there are still 
large uncertainties on the estimates presented in this pa- 
per, mostly due to the lack of precise knowledge on the 
fragmentation functions. Nevertheless, those uncertain- 
ties could be considerably reduced by using the data on 
hadron production collected at RHIC running with unpo- 
larized beams. Including these data on a global fit allows 
to define a safe procedure to extract the fragmentation 
functions needed to compute the measurable asymme- 
tries in the polarized case. 
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FIG. 8: NLO asymmetries forir + production at PHENIX with 
yfs = 200 GeV. 
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FIG. 9: The same as in FiglA but for 7r~ production. 



V. CONCLUSIONS 



We have computed the next-to-leading order QCD cor- 
rections to the spin-dependent cross section for single- 
inclusive hadron production in hadronic collisions. The 
calculation is implemented in a Monte-Carlo like code 
that allows to compute any infrared-safe quantity corre- 
sponding to one-hadron production to NLO accuracy, for 
both polarized and unpolarized collisions. 

Using the code, we have investigated the phenomeno- 
logical implications for pion production at RHIC. We 
find that the scale dependence is considerably reduced 
when the NLO corrections are included. Also the correc- 
tions are found to be non-trivial: K factors are larger for 
the unpolarized cross section than for the polarized one, 
resulting in a reduction of the asymmetry at NLO. The 
possibility of looking at charged pions in the final state 
is also studied, finding that 7r + and ir° are the most 
convenient channels for the determination of the polar- 
ized gluon density at values of x between 0.05 and 0.3. 
One hadron production at RHIC provides an excellent 
tool for the extraction of Ag, being possible to obtain 
a reasonable estimate of the size of the polarized gluon 
distribution just from the first run with longitudinally 
polarized beams to be performed in the next few months. 



Note Added 

As we completed this manuscript, we became aware of 
a calculation for the same process being performed by 
B.Jager, M.Stratmann and W.Vogelsang |53|. We have 
compared our results numerically finding good agree- 
ment. 
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